Vortex Dynamics in an Annular Bose-Einstein Condensate 
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We theoretically show that the topology of a non-simply-connected annular atomic Bose-Einstein 
condensate enforces the inner surface waves to be always excited with outer surface excitations and 
that the inner surface modes are associated with induced vortex dipoles unlike the surface waves 
of a simply-connected one with vortex monopoles. Consequently, under stirring to drive an inner 
surface wave, a peculiar population oscillation between the inner and outer surface is generated 
regardless of annulus thickness. Moreover, a new vortex nucleation process by stirring is observed 
that can merge the inner vortex dipoles and outer vortex into a single vortex inside the annulus. 
The energy spectrum for a rotating annular condensate with a vortex at the center also reveals the 
distinct connection of the Tkachenko modes of a vortex lattice to its inner surface excitations. 
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An atomic Bose-Einstein condensate (BEC) has been 
considered as an excellent many-body macroscopic quan- 
tum system with experimentally controllable parameters 
such as number density, inter-particle interaction, spin, 
confining potential and so on. Recently, topology of BEC 
also becomes to be important as in many different fields 
of physics. For example, the realization of a circular 
waveguide for an atomic BEC [1, 2] was followed by 
the achievement of a topologically non-simply-connected 
atomic condensate in a magnetic toroidal trap with a plug 
laser. It shows the first observation of atomic superfluid 
persistent current [3[ which was further refined to stay 
longer than 40 seconds in an all optical trap [4(. 

In torus BECs, there have been many interesting phe- 
nomena related with topological aspect: vortex lattice 
formation [EL @ , dynamics of vortices and solitons under 
stirring (7H1QI , and dynamics of vortex dipoles [Iol-[l2j . In 
order to understand the formation and decay of vortices 
in a superfluid related with the toroidal persistent cur- 
rent, the knowledge of low- lying energy spectrum and dy- 
namics is essential. When the torous BEC has flat shape 
with strong confinement along one direction, a toroidal 
condensate has two distinguishable surfaces, inner and 
outer surfaces (called as an annular BEC hereafter). The 
normal mode spectrum calculations for such an annular 
BEC available currently, however, have been done only 
within quasi- ID regime where the azimuthal degree of 
freedom is meaningful fl3j |. while in the experiments ra- 
dial degree of freedom cannot be ignored [3|, |4| . 

In this Letter, we calculate numerically the normal 
mode dynamics of an annular BEC within mean field 
regime considering the azimuthal and radial degrees of 
freedom fully compatible with current experiments d, Q . 
We have found that the excitations on the inner sur- 
face of a non-rotating annular BEC are associated with 
vortex dipoles rather than monopoles and that such in- 



ner surface excitations cannot be excited by themselves. 
Those aspects suggest topologically distinct phenomena 
of annular BECs when they are stirred, e.g., population 
oscillations between inner and outer surfaces and topo- 
logical excitations by vortex dipoles. A rotating annular 
BEC with nontrivial angular momentum is achieved by 
imposing quantized vorticity that penetrates the hollow 
region of the annulus. It is also found that such penetrat- 
ing vortices, which would give vortex lattice excitations, 
Tkachenko modes, if they were in a simply-connected 
BEC, play a crucial role in determining the inner sur- 
face excitations of a rotating annular condensate. 

To describe the normal mode excitations, we solve 
Bogoliubov-de Gennes equations [3, EH numerically us- 
ing finite element technique with Hermite polynomial ba- 
sis [18j ; Cuj—g^Vj = foujjUj and Cvj — g^Q 2 u 
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coupling constant is g 
ing length a, /i the chemical potential and the ground 
state wave function of the condensate, and ujj 

represent the wave functions and eigenenergy for the j-th 
quasiparticle excitation, respectively, and /i are de- 
termined by solving time-independent Gross-Pitaevskii 

equation (GPE), (-4^ + ^tr(r) + <#o| 2 - /i) = 
where is normalized by the total number of particles, 
N = fdr |^o | 2 - Strong confinement along the z-axis, 
uo xy <C uj z , is assumed. All values here are based on the 
dimension of the trap potential [16]. Fig.[T](b) is the den- 
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FIG. 1: (Color online) (a) Torus trapping potential, (b) The 
radial density distribution of BEC for various torus radii. 



FIG. 2: (Color online) Angular momentum (L z ) vs. energy 
(cj) curves for different ro values for a non-rotating conden- 
sate: r = (a), 2 (b), 4 (c), 6 (d) and 7 (e). 



sity of condensate along the radial direction for N = 5000 
and g = 0.1 for different ro values. 

Normal mode dynamics can be visualized by noting 
that ipj(t) = ipo + yf n^ (uje~ lUJj — v*e tuJj ) satisfies time- 
dependent GPE LL5|. Fig. 2 shows the energy vs. 
angular momentum dispersion relations of a BEC for dif- 
ferent ro values: (a) ro = and (b) ro = 2 correspond 
to simply-connected BECs, (c) ro = 4 for the transition 
region, and (d) and (e) (ro > 4) are for annulus cases. 
We categorize the normal modes using a pair of the ra- 
dial (n) and angular (m) quantum numbers, i.e., (n,m) 
mode, as in a simply-connected case; n and m represent 
the number of nodes along the radial and angular direc- 
tion respectively [17]. Here we notice that the energy 
dispersion curves as functions of m for either simply con- 
nected or annular condensate show well separated classes 
of normal modes with quantum number n's. 

First, we have found that the surface waves with the 
lowest radial quantum number (n = 0) cannot be gener- 
ated along the inner surface of an annulus even though an 
additional surface forms on its inner side. As ro changes 
from to 7 crossing the critical point at ro ~ 4, every nor- 
mal mode can adiabatically be transformed. The modes 
in the lowest curve (n = 0) in Figs. [2] (d) or (e) correspond 
to normal modes with longitudinal density wave propa- 
gating along the radial direction. Especially with higher 
|m|, n = density waves appear mainly as surface waves 
mostly on the outer surface as in the simply-connected 
case [17]; Fig. [3] (a) shows the contour plot of a snap- 
shot for the density of (0, 14) normal mode. Fig. [3] (d) 
represents the inner and outer edges if the annulus were 
to be unwound with the corresponding induced vortex 
singularity positions. Another schematic plot, Fig.[3](g), 
shows two density plots along the radial direction at two 
different angular positions A and B in Fig. [3] (a) which 
shows the density fluctuation at a certain fixed angular 
position. 

Next, we show that, unlike the case for the lowest radial 
mode (n = 0), both inner and outer surfaces host surface 
wave excitations for higher radial modes. In Fig. El (d) or 
(e), the next lowest curve (n = 1) corresponds to out-of- 
phase surface excitations with both inner and outer sur- 
face waves, of which wave nodes are placed alternatively 
along the angular direction. The corresponding modes in 
a simply-connected case, Fig. El (a) or (b), are those with 
one radial node including breathing mode [17j. The (1, 0) 
breathing mode in a simply-connected case has the well- 
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FIG. 3: (Color online) Contours for snapshots of \ipj(t)\ 2 
for modes with (n, m) = (a) (0,14) (b) (1,3) and (c) (2,3). 
Darker blue contour lines represent higher densities while red 
lines are drawn to clarify the density edges. Red (blue) dots 
represent induced (anti)vortices associated with the excited 
normal modes propagating conterclockwise and clarify where 
the wave nodes are. Contours (a), (b), and (c) are schemat- 
ically unwound in (d), (e) and (f) with vortex singularities 
propagting leftwards, respectively. Cross sectional density 
fluctuations along the lines A and B in (a), (b), and (c) are 
schematically drawn in (g), (h) and (i), respectively. 

known many-body breathing frequency 2uj xy [19]. One 
can see that such a breathing mode would transform 
into a radial sloshing mode centered at r = ro as the 
condensate gets annular-shaped and hence the frequency 
becomes corresponding radial trap frequency, uj xy . 

Figure [3] (b) shows the snapshot of (1,3) mode. It 
should be emphasized that the inner surface wave is ac- 
companied by vortex- ant ivortex pairs, vortex dipoles, dif- 
ferently from the outer surface wave that is accompanied 
by induced vortices or antivortices depending on the di- 
rection of the wave. The critical difference between the 
two surface waves is that the associated vortex singulari- 
ties for an inner surface wave are enclosed by the annulus 
of condensate while those for an outer one are not. Due 
to this, the induced vorticities on the inner surface will in- 
evitably generate supercurrent, a topological excitation, 
while those on the outer surface do not. By introduc- 
ing a vortex dipole inside the annular BEC the vortex or 
antivortex whichever is closer to the inner surface makes 
a propagating node in the surface wave while the other 
can annihilate the generated supercurrent. This prop- 
erty, which is a key finding in this work, is attributed 
to the topological aspect of annular geometry and the 
vortex excitations. A schematic plot, Fig. [3] (e), shows 
that an out-of-phase surface wave is in fact a transverse 
density wave. Furthermore, Fig. [3] (h) clarifies that the 
number of radial nodes (n) should be one for this specific 
radial center-of-mass dynamics, i.e., radial sloshing. 

The third curve (n = 2) in Fig. El (d) or (e) corre- 
sponds to in-phase excitations with both inner and outer 
surface waves, of which wave nodes are placed at the 
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same angular positions. The corresponding modes in a 
simply-connected case, Fig. [2] (a), are modes with two ra- 
dial nodes pi). As ro increases, modes with n = 2 make 
radial breathing dynamics centered at r — ro so that the 
frequency converges to near 2uo xy for m = 0. The (2,3) 
mode [Fig. [3] (c)] also shows vortex dipoles for the in- 
ner surface excitation. Figure [3] (f) is a schematic plot 
showing more clearly that an in-phase surface wave is a 
longitudinal density wave while Fig. [3] (i) shows that the 
number of radial nodes is two. In general, modes with 
higher n show similar surface waves, in-phase for even n's 
and out-of-phase for odd n's while they have additional 
radial nodes. 

Having categorized all the normal modes, it is noted 
that an inner surface wave should always be excited with 
an outer surface wave, either in-phase or out-of-phase. 
This implies that, if one drives an inner surface wave by 
stirring, it will not stay steady. It is rather a combina- 
tion of different in-phase and out-of-phase states which 
will generate beating, hence showing population oscilla- 
tion with outer-surface wave independent of the annu- 
lus thickness. In order to confirm this, we use a mod- 
ified trap potential mimicking inner surface mode such 
as V e = Vtr (1 + ecos(k6)) when r < ro and otherwise 
V e = Vt r (k is a non-zero integer). We obtain the ground 
state ^ e for V e that has an inner surface density modula- 
tion. Then let it evolve in time in Vt r after turning off e 
to be zero. Our simulation for an annular BEC of ro = 7 
shows that the inner surface density modulation propa- 
gates through the annulus reaching the outer surface and 
then back onto the inner, making a population oscillation 




FIG. 4: (Color online) The density fluctuation (n(t) = 
^(t) 2 — ^o [17]) of the inner (red) and outer (black) surface at 
— as a function of time when either only (a) the inner or 
(b) outer is excited at t = 0. Thin blue lines in (a) are guides 
for eyes. Snapshot of ip{t) — -y/pe 1 ^ under gaussian stirrer 
(white dot) with slow (c) and fast (d) speed. The height 
of the surface represents ^fp while the color <j). Red (black) 
circles represent ( ant i) vortices around which (j) has ±27r wind- 
ing. The black arrows in (c) and (d) denote the motions of 
vortices. 



between the two surface waves [Fig. |4] (a)]. On the other 
hand, we have also confirmed that an initial density mod- 
ulation on the outer surface makes a density fluctuation 
only there without any penetration at all [Fig. H(b)]. 

Another peculiar topological feature is in the nucle- 
ation process of vortices that can introduce a net topo- 
logical winding number within the annulus. As shown in 
Fig.[3^c), the annulus has a vortex- ant ivort ex pair inside 
while a vortex outside along the same angular position. 
If the annulus is stirred for inner and outer excitations 
to merge, they will annihilate themselves leaving a vor- 
tex inside the annulus, a topological excitation. Such 
vortex nucleation procedure cannot occur in a simply- 
connected case where phase singularities are typically 
generated outside and then get into the condensate via 
surface excitations. In order to confirm this senario, we 
have simulated stirring an annular BEC (ro = 10) with 
a narrow gaussian stirrer on the inner surface. Initially, 
the ground state is obtained with a gaussian stirrer fixed 
on the inner edge and then the speed of the stirrer is 
ramped up along the inner surface adiabatically at the 
angular acceleration of 10~ 5 uj^y / tt 2 . Figure 4 (c) shows 
the emergence of the first vortex dipole inside the annulus 
followed by the decay of the antivortex part of the dipole 
through the density dip near the stirrer leaving a vortex 
within the annulus generating supercurrent. Figure 4 (d) 
shows a similar case with two vortices already within the 
annulus at higher stirrer speed; a vortex dipole inside the 
annulus emerges together with a vortex outside followed 
by pair annihilation of the antivortex part of the dipole 
and the vortex outside on the rim of the annulus adding 
another vortex within the annulus. 

Now, let us look into the dynamics of an annular BEC 
with perpetual supercurrent; how the low- lying dynam- 
ics is affected by the current-carrying vortex at the torus 
center. Our calculations have been done in the lab frame. 
Figure [5] is the rotating-case counterpart of Fig. [2] with 
a single vortex at the center penetrating the BEC. Fig- 
ure [5] (a), a simply-connected case, is similar to Fig. [2] (a) 
except that it has a negative energy mode at m = — 1. 
This mode corresponds to the typical vortex precession 
that is the building block for Tkachenko modes in the 
case of vortex lattices [20]. It is interesting to note that 
Fig. [5] (e), an annulus case when compared to Fig. [2] (e), 
does not show such additional mode even though one 
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FIG. 5: (Color online) L z vs. uj for different ro values for a 
condensate with a vortex at the center: ro = (a), 2 (b), 4 
(c), 6 (d) and 7 (e). Red rectangle corresponds to a vortex 
precession mode at m — —1. 
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FIG. 6: (Color online) (a) Temporally successive (from left to 
right) snapshots of (0, —1) mode for a rotating simply-connect 
condensate and (b) those of (1,-1) mode for a annular con- 
densate with perpetual current. Blue dashed circles show the 
positions of ground state BEC. The red dot for each plot is 
the vortex that gives perpetual current to the BEC while the 
blue dot represents the induced antivortex accompanied by 
the m = — 1 excitation. 



can still make adiabatic deformation from Fig. [5] (a) to 
(e). We confirm that during the deformation, the vor- 
tex precession mode with a negative energy in Fig. [5] (a) 
adiabatically shifts to (0,-1) mode in Fig. [5] (e). 

This unusual adiabatic stabilization of a vortex preces- 
sion mode is made possible by the changes of the topol- 
ogy. As ro increases deforming the simply-connected con- 
densate into an annular one, the negative energy quasi- 
particle mode that makes the vortex state energetically 
unstable as a ground state acquires positive energy at 
ro ~ 0.8 because of the vortex pinning by the central po- 
tential bump. Once density hole is formed in the middle 
of the BEC with a large enough ro the precessing vor- 
tex behaves like the induced vortex singularities for an 
outer surface wave exciting the corresponding inner sur- 
face wave. This eventually forms annular density wave 
yielding (0,-1) mode. 

On the other hand, a mode with (0,-1) in a simply- 
connected condensate shows center-of-mass motion that 
makes a clockwise circular sloshing dynamics around a 
small circle with the vortex sympathetically following the 
density maximum [Fig. [6] (a)]. Such a sympathetic mo- 



tion of the vortex is a resonant behavior against the col- 
lective hydrodynamic excitation [17j. As the condensate 
transforms into an annular shape, this vortex motion also 
changes into a inner surface wave such that the nodes of 
inner and outer surface wave are out of phase, which cor- 
responds to (1,-1) mode [See Fig. 0(b)]. It is important 
to note that only m = — 1 modes can make resonant dy- 
namics with the singly-quantized vortex at the center of 
a rotating simply-connected BEC. With such resonant 
dynamics, (n, — 1) mode in a simply-connected BEC adi- 
abatically deforms into (n + 1, — 1) mode in an annular 
BEC. This resonant behaviour lowers the energy of the 
normal mode so that the m = — 1 mode has the lowest 
frequency for each family of n in a rotating annular BEC 
with a singly-quantized vortex at the center. This can be 
generalized to a rotating annular BEC with a vortex with 
higher winding number m v > 1 such that a mode with 
(0, — m v ) shows the characteristics of hydrodynamic sur- 
face excitations on the inner surface with m v nodes while 
it can be adiabatically transformed from the Tkachenko 
mode of a vortex lattice in a simply-connected case. 

In summary, our study on the Bogoiubov-de Gennes 
excitations of annular BEC reveals that the additional in- 
ner surface excitations are analytical continuations of col- 
lective excitations of a simply-connected condensate and 
that the topological aspect of the annulus forces inner 
surface waves associated with vortex dipoles. These as- 
pects give rise to interesting topological phenomena such 
as population oscillation and new vortex nucleation pro- 
cess. It is further found that, with perpetual current, the 
inner surface excitations have distinct connection to the 
the Tkachenko modes in a simply-connected condensate. 

Note added. Recently, we learned of a related 



preprint [21]. 
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